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ABSTRACT 

In a search for common proper motion companions using the VISTA Hemisphere Survey (VHS) and 
the 2MASS catalogs we have identified a very red (J — Kg = 2.47 mag) late-L dwarf companion of 
a previously unrecognized M dwarf VHS J125601.92-125723.9 (hereafter VHS 1256-1257), located at 
a projected angular separation of 8706±(y.'03. In this work we present a suite of astrometric, pho¬ 
tometric, and spectroscopic observations of this new pair in an effort to confirm the companionship 
and characterize the components. From low-resolution (i?^ 130-600) optical and near-infrared spec¬ 
troscopy we classified the primary and the companion as an M7.5 ± 0.5 and L7 ± 1.5, respectively. 
The primary shows slightly weaker alkali lines than field dwarfs of similar spectral type, but still 
consistent with either a high-gravity dwarf or a younger object of hundreds of millions of years. The 
secondary shows spectral features characteristic for low surface gravity objects at ages below several 
hundred million years, like the peaked triangular shape of the i7-band continuum and alkali lines 
weaker than in field dwarfs of the same spectral type. The absence of lithium in the atmosphere of the 
primary and the likely kinematic membership to the Local Association allowed us to constrain the age 
of the system to the range of 150-300 Myr. We report a measurement of the trigonometric parallax 
71 = 78.8 ±6.4 mas, which translates into a distance of 12.7 ±1.0 pc; the pair thus has a projected 
physical separation of 102 ± 9 AU. We derived the bolometric luminosities of the components and 
compared them with theoretical evolutionary models to estimate the masses and effective tempera¬ 
tures. For the primary, we determined a luminosity of log(Lboi/A 0 ) = —3.14 ± 0.10, and inferred a 
mass of 73^15 Afjup at the boundary between stars and brown dwarfs and an effective temperature of 
2620± 140 K. For the companion we obtained a luminosity of log(Lboi/± 0 ) = —5.05±0.22 and a mass 
of 11.2l®'g Mjup placing it near the deuterium-burning mass limit. The effective temperature derived 
from evolutionary models is 880^)^20 about 400-700 K cooler than the temperature expected for 
field late-L dwarfs. 

Subject headings: stars: brown dwarfs - stars: imaging - infrared: planetary systems - stars: individual 
(VHS J125601.92-125723.9) 


1. INTRODUCTION 

Very low-mass stars and brown dwarfs directly imaged 
around stars constitute an important group for studies of 
the properties of substellar objects. We can infer the dis¬ 
tance and metallicity of the companion from the brighter, 
more easily characterizable primary star. Most impor¬ 
tantly, we can constrain the age of the system and thus 
overcome the intrinsic degeneracy between mass and age 
for the temperature and luminosity of objects below the 
hydrogen burning limit. This allows for a more thorough 
characterization and offers the opportunity to better un¬ 
derstand their physical properties a nd to test evolution- 
ary and atmospheric models (e.g., Pinfield et al. 2006 
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Faherty et al.|2010 1 . There are about 1500 objects spec- 


traliy classihed as M7 and later, the vast majority is 
found to be isolated, single sources that appear to have 
spectral energy distributions, photometric colors, and 
kinematics consistent with a field population with ages in 
the r ange from 2 to 8 Gyr pTnapp et al.|2004 Cruz et al. 


2007 


Faherty et al. 2009|)' About 120 uitracool dwarfs 


(late-ivl, L, 1’, and Y dwarfs) have been confirmed as 


components of binary or multiple systems (Faherty et al. 
2010 Deacon et al.|2014). These objects discovered over 


the last two decades have exhibited a large diversity in 
ages, atmospheric properties, and chemical compositions. 
Surveys aimed at identifying low-mass companions of 


young stars (e.g., | Wahhaj et al.|20Tl| [Bowler et al.|2013 


vounff 

star clusters (e.g.,|Beiar et al. 120111 

Pena Ramirez 

et al. 

2012 

Zapatero Osorio et al. 

20T4 

M, as well as 

searches tor 

L and 1' dwarfs in the 

ieid using large sky 


area surveys (Kirkpatrick et al. 2006 Cruz et al. 20091 
have revealed a number of ultracool dwarfs that exhibit 
photometric and spectral properties different from typ¬ 
ical late-type objects of the field population. Some of 
these peculiarities have been attributed to low surface 
gravities and cloudy atmospheres, occurring at the early 
stages of evolution, below several hundred million years 
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Fig. 1. — VISTA YJHKs band images of the new common proper motion pair VHS 1256-1257. The identified companion is located 
at a position angle of 218?ldi0?2, with a projected angular separation of 8706di0703 which corresponds to ~ 102 AU at the determined 
parallactic distance of 12.7 pc. The field of view is 30''x30'', with north up and east to the left. 


(ICruz et al. 

2009 

Allers & Liu 

2013 

Zapatero Osorio 

et al. 2014a 

. Young L dwarts, either tree-hoatmg or as 


companions to stars, were found to share similar charac- 
teristics. They have very red colors {J — Ks>2 mag), 
J-band absolute magnitudes fainter than their old, field 
counterparts, and show distinctive spectral features such 
as, for example, a sharply peaked, triangal-shaped con¬ 
tinuum in the H band and weaker sodium and potassium 
lines. 


Recent studies have revealed a strong resemblance be¬ 
tween the young L dwarfs and directly imaged plane- 


tary mass companions (e.g., 2MASS 1207-39 b, 

Chau- 

|vin et al. 

20051 HR 8799 bcde, 

Marois et al.||2008 

■pm 

CJ 504 b 

Kuzuhara et ai.||2013 

S Pic b, Lagrange et al.| 

2009). T 

ley have similar near-infrared (near-IR) colors 


and absolute magnitudes, overlapping effective tempera¬ 
ture regimes of ~ 1000-150 0 K and masses of a few to a 


few tens of Jupiter masses 

(Barman et al. 2011a Bowler 

et al. 

2013 Faherty et al. 

2m^' 

2014 Liu et al. 

20n). 


bpectroscopic and photonietric studies ot young substel- 
lar objects can provide information on the physical prop¬ 
erties of gas giant exoplanets found by transit and radial 
velocity surveys, in particular, on the characteristics and 
composition of their complex atmospheres. 

In this work we present the identification and charac¬ 
terization of a nearby, young binary system with com¬ 
ponents that we classify as M7.5 ± 0.5 and L7 ± 1.5. In 
Section 2 we describe the data and search method that 
led to the identification of the pair. Section 3 contains 
the description of observations aimed at confirming their 
companionship and characterization. In Section 4 we dis¬ 
cuss their physical properties. We classify their spectral 
types, and determine the distance and radial and space 
velocities. We constrain the possible age of the system 
and estimate the luminosities, masses, and effective tem¬ 
peratures of the two components. Final remarks and 
conclusions are presented in Section 5. 


2. IDENTIFICATION OF THE SYSTEM 
2.1. VISTA Hemisphere Survey Data 

The reported common proper motion pair was identi¬ 
fied using th e VISTA Hemisphere Survey (VHS) data and 
the 2MASS (jSkrutskie et al.l2006|. The VHS is a near-IR 
ESO public survey designed to map the entire southern 
hemisphere in the J and broad-band filters with av- 
erage 5a depths of J = 20.2 mag and Kg = 18.1 mag 
(McMahon et al. 2013). In some particular areas also 
observations in the Y and H bands are also performed. 


The 4-m VISTA telescope (Emerson|2001 Emerson et al. 
2004) has operated since 2009 at ESO's Oerro Raranal 
O bservatory in Chile and has thus far covered about 8000 
deg^ of the sky. It is equipped with a w ide-field infrared 
camera VIRCAM (|Dalton et al. 20061 composed of 16 
Raytheon detectors 2048x2048 pixel array each, with a 
mean plate scale of 0'.'34, giving a field of view of 1.65° 
in diameter. 

The VHS images are processed and calibrated auto¬ 
matically by a dedicated science pipeline implemented by 
the Cambridge Astronomical Survey Unit (CASH). Stan¬ 
dard reduction and processing steps include dark and sky 
subtraction, flat-field correction, linearity correction, de¬ 
stripe, and jitter stacking. For a detailed description we 
refer the reader to the CASH webpage http://casu. 
ast■cam.ac.uk/s urve ys-projects/vist a as well as to 
Irwin et al. (2004) and Lewis et al. (2010). 

'I'he photometry provided in the VHS catalog is cal¬ 
ibrated using the 2MASS magnitudes of all matching 
stars converted onto the VISTA system applying color 
equationt]^ which include terms accounting for the in¬ 
terstellar reddening. Photometric calibrations are de¬ 
termined to an accuracy of l%-2%. The astrometric 
solution for VHS observations is obtained through the 
CASH pipeline, using the 2MASS point source catalog. 
The objects on the catalogs extracted from each detec¬ 
tor are matched to their counterparts in 2MASS using 
a correlation radius of 1". Because 2MASS has a high 
degree of internal consistency it is possible to calibrate 
the world coordinate system of VISTA images to relative 
accuracy better than O'.'l. 


2.2. Search Method 

Using VHS data, we carried out a search for high 
proper motion objects, by cross-matchi ng VHS sources 
with the 2MASS Point Source Catalog (Skrutskie et al. 
2006). The search focused on objects that had moved at 
least 2"and a maximum of 30"from 2MASS to VHS. The 
time baseline between the two surveys is typically about 
12 yr, which gives proper motions of approximately 0.15- 
3.0 arcsecyr”^. Over the common area of ~5000 deg^ 
between 2MASS and VHS we found more than 6000 
objects with J magnitudes in the range ll-17mag and 
proper motion higher than 150 mas yr~^. We have cross - 
correlated this catalog with WISE ) Wright et al.||2010 ) 


so that near and mid-infrared inforniation is available tor 


® http://casu.ast.cam.ac.uk/surveys-projects/vista/ 

technical/photometric-proparties 
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each target. Most of the high proper motion objects are 
relatively nearby M dwarfs with estimated photometric 
distances <50 pc (A. Perez-Garrido et ah, in prepara¬ 
tion) . We have searched for common proper motion pairs 
and multiples among these objects with proper motion 
consistent within 40 masyr“^ in both right ascension and 
declination (/r„, ns)- VHS 1256-1257 system was one of 
the identified candidates. The pair was observed with 
the VISTA YJHKs filters and with a seeing of I'.'l on 
2011 July 1. The primary and the secondary were cat¬ 
aloged in 2MASS with designations 2MASS J125602.15- 
125721.7 and 2MASS J125601.83-125727.6, respectively. 
The primary was also listed as SIPS 1256-125 7 in the 


sample of low-ma ss stars with p, > Ol'l yr ^ from Deacon 


& Hambly (20071, with measured proper motion ampli¬ 


tude fi = u! 357 yr“^ and the position angle of the proper 
motion vector of 244.76°. 

VISTA YJHKs images of the pair are presented in 
Fig.[T] and the photometry is given in Table The com¬ 
panion is located at a projected angular separation of 
8'.'06±0'.'03, at a position angle of 218?1±0?2. This sep¬ 
aration corresponds to a projected orbital separation of 
102 ± 9 AU at the estimated distance of the system (Sec¬ 
tion 4.2). The two components share a common proper 
motion, which significantly differs from the proper mo¬ 
tion of background stars as shown in Fig.[^ The /i^ cos 6 
and fis measured from the VHS and 2MASS positions of 
the sources, were —270 ± 17, —185 ± 14 and —292 ± 27 
and —212 ± 23 masyr“^, for the primary and secondary, 
respectively. The uncertainties in proper motion of each 
component correspond to the rms of /io, cos S and fis of 
the background stars within a one-degree radius around 
the primary, and magnitudes similar to that of the given 
component. Time baseline between the two epochs is 
12.3 yr. The primary has J = 11.02 ± 0.02 mag and 
J — Ks = 0.97 mag (2MASS photometry, VIRCAM is 
out of the linear range) The secondary is roughly six 
magnitudes fainter in the J band (J = 17.14±0.02 mag) 
and has a very red J — Kg color of 2.47 mag indicating 
its significantly cooler type and lower mass with respect 
to the primary. 


3. FOLLOW-UP OBSERVATIONS AND DATA 
REDUCTION 

3.1. NTT/SofI Near-infrared Spectroscopy and 
Imaging 

To measure the infrared spectral types of the 
VHS 1256-1257 components we performed follow-up 
near-IR spectroscopy using the Son of ISAAC (Sofl) 
spectro-imager installed on the 3.6 m New Technology 
Telescope (NTT) on 2014 March 12. Sofl is equipped 
with a Hawaii HgCdTe detector with 1024x1024 pixel 
array offering a field of view of 4.9x4.9 arcmin^ with a 
Cf.'288 pixel scale. We have used blue and red grisms 
covering the 950-1640 and 1530-2520 nm range, com¬ 
bined with a slit of 1", orientated along the position 
angle of the components of the system. This instru¬ 
mental configuration provides a nominal dispersion of 
6.96 and 10.22 Apix“^ for the Blue and Red grating, 
respectively. To subtract the sky background we used an 
ABBA nodding pattern with an offset of 20"between the 
two positions. Individual exposure times were 600 and 
900 s for the red and blue grism, respectively. Right af¬ 



Fig. 2. — Proper motion diagram for the VHS 1256-1257 system. 
All correlated objects within one degree around the primary with 
J <17.5 mag are plotted with black dots. Components of the com¬ 
mon proper motion pair are labeled. Error bars correspond to the 
astrometric RMS of stars near the magnitudes of the components. 


ter VHS 1256-1257 we observed a telluric standard star 


van 


HD 112304 (V = 6.19 mag, AOV, Hpg et al. 2000. . 
Leeuwen]|2007 1. The sky conditions during the observa- 


tions were clear with sub-arcsecond seeing. 

The two-dimensional (2D) images were reduced and 
combined using the ESO Sofl pipeline run within the 
Gasgano environment. The traces of the primary and 
the companion spectra on the 2D images were clearly sep¬ 
arated and there was no contamination between the two 
sources. The two spectra of the VHS 1256-1257 system 
and of HD 112304 were then extracted using standard 
IRAE routines within the APALL package and wavelength 
calibrated via Xenon arc lines. The dispersion solution 


had an rms of 0.4 and 0.35 A for the blue and red part 
of the spectrum, respectively. The final resolution of the 
spectrum was 24 A (i? ^ 600). Telluric absorption lines 
were corrected, dividing the spectra by the AOV telluric 
standard HD 112304 observed at a similar airmass, and 
multiplying by a blackbody of a corresponding effective 
temperature of 9480 K. The obtained NTT Sofl spectra 
of VHS 1256-1257 primary and companion normalized at 
1.6/rm are displayed in Eigurel^ 

On 2014 April 27 we used Sofl in imaging mode, to 
perform J- and H-band observations at an additional 
epoch. We used a nine-position dither pattern with in¬ 
dividual exposures of 5 s. The weather was clear, with a 
seeing of 1'.'15-1'.'30. Raw images were reduced using the 
ESO Sofl Pipeline within Gasgano, which includes bias 
subtraction, flat-field correction, plane re-centering, and 
co-addition of individual frames to a final image. Images 
were employed for determining the trigonometric paral¬ 
lax of the system and for performing relative astrometry 
between the two components as described in Section 4.2. 
Table 1^ provides these measurements. 


3.2. GTC/OSIRIS Optical Spectroscopy 

We obtained low-resolution optical spectroscopy of 
VHS 1256-1257 with the Optical System for Imag- 
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TABLE 1 

Observation log of VHS 1256-1257. 


Obs. Date 
(UT) 

Tel/Instrument 

Mode 

Wavelength 

(/rm) 

Exp. Time 

(s) 

Seeing 

(") 

Slit 

(") 

Grating 

Scale 

("/Pix) 

Res. 

Power 

2011 Jul 01^ 

VISTA/VIRCAM 

Img 

YJHKs 

15, 15, 7.5, 7.5 

1.0 

- 

- 

0.34 

- 

2014 Mar 02*^ 

GTC/OSIRIS 

Spec 

0.48-1.00 

120 

0.57 

1.5 

R500R 

0.25 

320 

2014 Mar 12^ 

NTT/SoH 

Spec 

0.95-1.64 

4x600 

0.8 

1.0 

Blue 

0.29 

600 

2014 Mar 12^ 

NTT/SoH 

Spec 

1.53-2.52 

4x900 

0.8 

1.0 

Red 

0.29 

600 

2014 Apr 22, May 04*= 

NOT/ALFOSC 

Spec 

0.63-0.68 

900, 4x900 

1 

00 

d 

1.3 

Grism#17 

0.19 

4700 

2014 Apr 27“ 

NTT/SoH 

Img 

JH 

9x2 

1.0-1.3 

- 

- 

0.29 

- 

2014 May 14*’ 

GTC/OSIRIS 

Spec 

0.56-0.77 

3x600 

0.56 

0.6 

R2500R 

0.25 

3000 

2014 Jun 03^ 

GTC/OSIRIS 

Spec 

0.48-1.00 

3x600 

0.63 

1.5 

R300R 

0.25 

130 

2014 Jul 01, 09, 10*’ 

VLT/UVES 

Spec 

0.37-0.49 

3x3300 

0.5-1.7 

1.0 

CD2, 4 

0.22 

40000 

2014 Jul 01, 09, 10*’ 

VLT/UVES 

Spec 

0.56-0.95 

3x3300 

0.5-1.7 

1.0 

CD2, 4 

0.16 

40000 

2014 Jul 15*’ 

IAC80/CAMELOT 

Img 

VI 

300 

1.1 

- 

- 

0.30 

- 

2014 Jul 17^ 

WHT/ACAM 

Img 

Sloan i, z 

9x30, 9x30 

0.85 

- 

- 

0.25 

- 

2014 Jul 18^ 

WHT/LIRIS 

Img 

j 

5x9x2 

0.7 

- 

- 

0.25 

- 

2014 Dec 19, 20*’ 

IAC80/CAMELOT 

Img 

I 

6x300 

1.5, 1.2 

- 

- 

0.30 

- 


^ Both components were observed 
^ Only the primary was observed 


ing and low-inter mediate Reso lution Integrated Spec¬ 
troscopy (OSIRIS; |Cepa||2010|) mounted on the 10.4- 
m Gran Telescopio de Ganarias (GTC) telescope in La 
Palma. OSIRIS is equipped with two 2048x4096 Mar¬ 
coni GGD42-82 detectors offering a field of view approx¬ 
imately 7x7 arcmin^ with an unbinned pixel scale of 
0'.'125. We observed the primary star using two differ¬ 
ent gratings: R500R, which allowed us to measure the 
general spectral energy distribution at a resolution of 
i?~320, covering the 0.48-1.00 ^m range, and R2500R 
providing a resolution of i?~3000 in the 0.56-0.77 ^m 
region, which we used to study the Ha and Lii lines. 
Observations with R500R and R2500R were acquired on 
2014 March 2 and May 14, respectively. On 2014 June 
3, both the primary and secondary were aligned on the 
l75-width slit and were observed simultaneously using 
the R300R grating. This instrumental setup provided 
optical spectra with a resolving power of 130, which 
is optimal to maximize the signal-to-noise ratio of the 
faint companion data. Observations were performed in 
service mode (Table as part of a GTC filler program 
(number GTC5-14A; PI Lodieu). The nights were dark 
with a clear sky and a seeing of ^0'.'6. Details of the 
used instrumental configurations and exposure times are 
provided in Table [T] The spectrop hotometric standards 
HiltnerOOO, GD 153 , and Ross 640 (|H0g et al.|1998||2000, 
Pancino et al. 2012[) were observed with the R500 R+Gr' 


R250UR, and R300R gratings, respectively, on the same 
night of the scientific observations. Observations of the 
standards with R300R and R500R were done with the 
same grating combined with the Sloan z filter to correct 
for second -order contamination beyond 9 200 A (see pro¬ 
cedure in Zapatero Osorio et al. 2014a). Bias frames, 
continuurn lamp fiat fields, and JN eon-|- Aenon arc lamps 
were observed by the observatory staff during the after¬ 
noon preceding the observations. 

We reduced the OSIRIS spectra using routines within 
IRAK. We subtracted the raw spectra by a median- 
combined bias and divided by a normalized continuum 
lamp flat field. From the optimally combined 2D images 
we extracted the spectra using the APALL routine and cal¬ 
ibrated in wavelength with the lines from the combined 
arc lamp. To improve the detection of the secondary in 
observations with R300R we inserted both components 


in the slit and used the trace of the primary as a reference 
to combine the individual exposures. Gorrection of the 
instrumental response was done using the corresponding 
spectrophotometric standards observed during the same 
nights. The obtained spectra are displayed in Figures 
|4]and|6l 

3.3. NOT/ALFOSC Optical Spectroscopy 


For a first estimate of the radial velocity of the pri¬ 
mary star, intermediate-resolution optical spectroscopy 
was acquired for VHS 1256-1257 A using the Andalucia 
Faint Object Spectrograph and Gamera (ALFOSC) in¬ 
strument of the 2.5 m Nordic Optical Telescope (NOT) 
operating at the Observatorio del Roque de los Mucha- 
chos. La Palma. The camera uses a E2V back illuminated 
GGD42-40 chip with 2kx2k pixels with a scale of 0.19 
arcsecpix“^, providing a field of view of 6.4x6.4 arcmin^. 
We observed the primary in two different epochs, on the 
nights of 2014 April 22 and May 5. For the instrument 
setup we chose Grism#17 and a 173 slit providing a wave¬ 
length range of 6330.4-6853.6 A and a nominal resolution 
of 0.255 A pix“^. On the first night one useful exposure of 
900 s was obtained under the presence of variable clouds. 
On the second night four exposures of 900 s were taken, 
under clear sky conditions. The seeing was similar on 
both nights (078-171), and hence the final resolution of 


the spectra was also very similar: 1.5 A (i? ^ 4700). Dur¬ 
ing the same nights using the same instrumental config¬ 
uration, we also observed the star GJ 388, which is a 
M4.5V with a precise radial yelocitv determinatio n of 
Vr = 12.453 ± 0.066 kms“^ by Ghubak et al. (2012). 

Raw data were reduced using routines within the iraf 
environment. Two dimension images were bias corrected 
and fiat fielded using continuum lamps normalized using 
the RESPONSE routine. Spectra were optimally extracted 
using the APALL routine and were wavelength calibrated 
using HeNe lamps and a cubic spline function fit of the 


order of three, providing an rms better than 0.06 A. Fi¬ 
nal spectra were corrected by the instrumental response 
using the spectrophotometric standard HZ44. 


3.4. VLT/UVES Optical Spectroscopy 

To measure a more precise radial velocity and to inves¬ 
tigate in detail the spectral signatures of youth of the pri- 
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mary star, like the Lii line at 670.82 nm we also obtained 
a high-resolution spectrum of VHS 1256-1257 A using 
the Ultraviolet and Visual Echelle Spectrograph (UVES) 
mounted on the Kueyen unit of the ESO Very Large 
Telescope. UVES is a two-arm crossdispersed echelle 
spectrograph covering the wavelength range 300-500 nm 
(blue) and 420-1100nm (red), with the possibility to use 
dichroic beam splitters. The instrument is equipped with 
a single chip in the blue arm and a mosaic of two chips in 
the red arm. The blue CCD is a 2Kx4K, 15 /im pixel size 
thinned EEV CCD-44. The red CCD mosaic is made of 
an EEV chip of the same type and the MIT/LL CCID-20 
chip. Each arm has two cross disperser gratings work¬ 
ing in first spectral order; the typical order of separation 
is 10". We have used the standard dichroic mode set¬ 
ting, which covers the 565-950 nm spectral range in the 
red arm, with the central wavelength at 760 nm, and a 
373-499 nm range in the blue arm, centered at 437 nm. 
This configuration, with 1 arcsec slit, provides spectral 
resolution of i? ^ 40 000 (25 mA/pix). Observations were 
performed in service mode under the ESO DDT program 
293.C-5014(A) on 2014 July 1, 9 and 10, with clear sky 
conditions. The total integration time was 9900 s divided 
into three exposures of 3300 s. 

The spectra were reduced, extracted, and calibrated 
using the ESO UVES pipelines under the ESO Recipe 
Flexible Execution Workbench environment {Reflex, 
Freudling et al.|[2013 ). The reduction steps of the work- 
how executed by the software include creation of master 
flat and bias frames, with the corresponding corrections 
of science data and detection of the order positions on 
the detector. Subsequently, a wavelength calibration so¬ 
lution is obtained from an input arc-lamp frames. We 
have used the ThAr arcs acquired in both blue and red 
arms, during the afternoon preceding the observations. 
Instrument characteristics description provides a note of 
caution in that in the spectral region above 700 nm the 
ThAr lamp has some very bright Argon lines that satu¬ 
rate the CCDs, and the heavily saturated lines may con¬ 
tain remnants in the following exposures. Here the cal¬ 
ibration frames were taken more than four hours before 
the night, i.e. more than the typical time after which the 
remnants vanish. Then, the spectrum from each order is 
extracted and merged. Last, the flux-calibration of the 
science spectrum is carried out, using the appropriate in¬ 
strument response curve, in the blue part obtained from 
a standard star observation, and from the instrument 
master response curve in the red part. We have corrected 
the telluric absorption f eatures in the spectrum using the 
ESO Moleefit software (Smette et al.||2015 Kausch et al. 
20151. The results obtained using the U VES data are 
described in Section 4.4, with the final spectra presented 
in Figures [T0| and [g 

3.5. IAC80/CAMELOT Imaging 

On 2014 July 15 and December 19 and 20 we performed 
imaging observations of the primary using CCD Camera 
CAMELOT of the 0.8 m IAC80 telescope at the Teide 
Observatory on Tenerife. We aimed to obtain optical 
photometry in the VI (Johnson-Bessell’s) filters. The 
images were also used in the parallax determination. The 


TABLE 2 

Measurements and determined physical parameters of 
VHS 1256-1257 system 


Astrometry 

Primary 

Companion 

R.A.(J2000)^ 

12*'56“01=.922 

12*'56™01±586 

Decl. (J2000)^ 

-12°57'23'.'990 

-12'=57'30'.'310 

2MASS ID 

J125602.15-125721.7 

J125601.83-125727.6 

Separation (arcsec)^ 

8.06 ±0.03 

Separation (AU) 

102 ±9 

Position angle (deg)^ 

218.1 ±0.2 

pQ, cos 5 (masyr“^) 

-281.5 ±5.3 

-275.4 ± 5.3 

fis (rnasyr”^) 

-205.5 ± 15.2 

-198.4 ± 15.2 

Parallax tt (mas) 

78.8 ±6.4 

Distance d (pc) 

12.7 ± 1.0 

Vr (kms“^) 

-1.4 ±4.5 


U (kms“^) 

-9.4 ±2.0 


V (kms“^) 

-16.4 ± 3.0 


W (kms“^) 

-8.5 ± 3.9 


Photometry 

V (mag) 

17.759 ±0.059 


I (mag) 

Sloan i (mag) 

13.318 ± 0.039 

22.494 ± 0.315 

Sloan 2 (mag) 


20.095 ±0.090 

2MASS J (mag) 

11.018 ±0.023 

16.662 ±0.287 

2MASS H (mag) 

10.473 ± 0.023 

15.595 ±0.209 

2MASS Ks (mag) 

10.044 ± 0.021 

14.568 ±0.121 

VHS Y (mag) 

<11.72 

18.558 ±0.051 

VHS J (mag) 

<11.36 

17.136 ±0.020 

VHS H (mag) 

<11.02*= 

15.777 ± 0.015 

VHS Ks (mag) 

<10.42*= 

14.665 ±0.010 

WISE W1 (mag) 

9.880 ±0.023 

13.6 ±0.5 

WISE W2 (mag) 

9.658 ±0.021 

12.8 ±0.5 

WISE W3 (mag) 

9.390 ± 0.044 

> 11.8 

WISE W4 (mag) 

8.334 ±0.410 

>8.65 

Spectral Classification 

Optical 

M7.0±0.5 

L8.0± 2.0 

Near-IR J 

M8.0±0.5 

L8.0± 1.0 

Near-IR K 

M8.0±0.5 

L5.0± 2.0 

H 2 O index 

M7.6±0.4 

>L4.0 

H 2 OD index 

<L0 

L5.8± 0.8 

H 2 O-I index 

M6.9± 1.1 


H 2 O -2 index 

M7.9±0.5 


Adopted spectral type 

M7.5±0.5 

L7.0± 1.5 

Near-IR gravity class 

INT-G 

VL-G 

Physical Properties 

Age (Myr) 

150- 

-300 

log(Lboi/L0) 

Mass (Mjup) 

Tefl (K) 

-3.14± 0.10 

70 + 20 

2620 ± 140 

-5.05 ±0.22 

11.2l?;^ 

880ll« 

log(5) 

5.05 ± 0.10 

4.25l°;?5 


^ Measured using VHS 
YJHKs photometry 


images, epoch (MJD) = 55743.067635 
out of VIRCAM linear range 


in the V and I filters. On the night of July 15, right after 
the scientific target, we observed two standard star fields 
at similar airmass, containing eight standard stars used 
for the photometry calibration. Weather conditions were 
clear/photometric with seeing of 1'.'1-1'.'3. 

The data were reduced automatically using a dedi¬ 
cated PyRAF script that includes a standard bias and 
flat-field corrections procedure. We carried out aperture 
photometry with the phot task within the iraf envi¬ 
ronment. Photometric calibration of the instrumental 


with 07304 pixels corresponding to a 10.4x10.4 arcmin^ 

dard stars ( 

Landolt 

1992 

2009 

1 from the two observed 

field of view. Individual exposures of 300 s were obtained 

tieids. i’he uncertainty ot 

the calibration is ~U.U4 mag. 
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Fig. 3.— Low-resolution optical (0.5-l.Op.m, 120-320) and near-infrared (1.0—2.4pm5 i?~600) spectra of the primary (top) and 

companion (bottom) obtained using GTC/OSIRIS and NTT/Sofl. The flux is normalized at 1.6/im. The most prominent molecular and 
atomic features are indicated. Unlike most of old population field L dwarfs, the spectrum of VHS 1256-1257 b peaks in the H band, which 
has a triangle-shaped continuum characteristic for young, dusty L dwarfs. 


The measured VI magnitudes of the primary and other 
photometric data of the primary and the companion are 
listed in Tabled 

3.6. WHT/ACAM and LIRIS Imaging 

To complement the photometric information at differ¬ 
ent wavelengths, we also acquired optical images of the 
pair using the Sloan i and 2 broad-band filters. The final 
z-band image was also employed to constrain the trigono¬ 
metric parallax (Section 4.2.2). Observations were done 
using the ACAM camera at the 4.2-m William Herschel 


Telescope (WHT) on 2014 July 17. ACAM uses a 2kx2k 
pixels EEV CCD detector with 0.25 arcsecpix”^ scale, 
providing a circular field of view with a 8.3 arcmin^ di¬ 
ameter. We used a nine-point dither pattern with indi¬ 
vidual exposures of 30 s. On the same observing run on 
July 18 we have obtained J-band observations using the 
Long-slit Intermediate Resolution Infrared Spectrograph 
(LlRlS) spectrograph and imager on the WHT, for an 
additional astrometry epoch of the binary used to con¬ 
strain the parallax measurement. LlRlS uses a Ikxlk 
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Fig. 4. — Low-resolution GTC/OSIRIS optical (left) and NTT/SofI near-infrared (right) spectra of the primary VHS 1256-1257 plotted 
with black line, compared with late-M spectral templates (with labels indicating their type). Comparison spectra were smoothed to match 
the resolution of our data. The sources and references of used templates are described in Section 4.1.1. Spectra were normalized at 0.9 /rm 
at optical and at 1.6 fira at near-IR and offset by a constant for display. 


HAWAII detector with a pixel scale of 0.25 arcsecpix”^, 
yielding a field of view of 4.27x4.27 arcmin^. We used a 
sequence of five dither patterns of 9 positions and 2 s in¬ 
dividual exposures, giving a total exposure time of 90s. 
Weather conditions were clear during the observations, 
with average seeing of 0'.'85 on the first and 0'.'66-0'.'68 
on the second night. 

We carried out the data reduction of the LIRIS im¬ 
ages with the LIRIS data reduction packag^ The pro¬ 
cedure consisted of subtracting the sky made from the 
dithered images, correction for fiat field, vertical gradient 
observed on the detector and the geometrical distortion. 
ACAM images were reduced with standard procedures 
using IRAF. Master sky flat frames were obtained from 
science images by combining by the median the individ¬ 
ual exposures of the dither sequence. Corrected images 
were aligned and average combined. The aperture pho¬ 
tometry of the companion was obtained using the phot 
task within iraf. To calibrate the instrumental mag¬ 
nitudes we measured the photometry of stars in a field 
observed after VHS 1256-1257, under similar weather and 
airmass conditions. We selected 10 point sources with a 
good signal-to-noise ratio and av ailable measurem ent in 
the SDSS Photometric Catalog (Ahn et al. 201^ The 
obtained i and 2 magnitudes are listed in 'i'abie|^ The 
primary was saturated in the individual exposures. 


3.7. WISE Data 

The Wide-field Infrared Survey Explorer [WISE] 
Wright et al. 2010) conducted a survey of the entire 
sky in the 3.4, 4.6, 12 and 22 yrm mid-infrared band- 
passes (hereafter Wl, W2, 1V3 and WA). The HTS'i? All- 
Sky Source Catalog (see explanatory supplement docu¬ 
ment by Cutri et al. http: //wise2 . ipac . caltech. edu/ 
docs/release/allsky/expsup/ includes the photome¬ 
try of VHS 1256-1257 in the four bands, but does not 
detect VHS 1256-1257 b. The kWWISE Source Catalog 
contains the deblended photometry of both the primary 
in the four bands and the secondary in the Wl and W2, 


^ WWW.ing.iac.es/Astronomy/instruments/liris/liris_ql. 
html 




X (nm) X (nm) 



Fig. 5.— Comparison of the Nal and Kl lines of the primary 
VHS 1256-1257 in the optical (top panels) and near-IR (bottom 
panel) with standard field M7—M8 counterparts and with young 
low-gravity M8-type dwarfs Teide 1 and 2M 1207 A. All spectra 
per panel were conveniently normalized at a certain wavelength, 
and lines were convolved to the same resolution for a proper com¬ 
parison. These absorption features are slightly less intense than in 
field objects of similar spectral type. 


but the later is contaminated by a diffraction spike of 
the bright star as indicated in the confusion flag. In or¬ 
der to obtain the correct values, we performed aperture 
photometry of VHS 1256-1257b using daophot routine 
in the Wl and lT2-band images after PSF subtraction 
of the primary using the scaled PSF of nearby stars in 
the field as references. A small aperture of 5 pixels was 
used and instrumental magnitudes was transformed into 
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apparent magnitudes using the photometry of VHS 1256- 
1257. The KWWISE Source Catalog photometry of the 
primary and the derived photometry of the secondary are 
given in Table 

4. PHYSICAL PROPERTIES OF THE SYSTEM 
4.1. Spectral Types and Spectral Characteristics 

In Figure we plot the merged low-resolution optical 
(0.5-I.0/rm) and near-IR (1.0-2.4p,m) spectra of each of 
the two components of VHS 1256-1257. The spectral en¬ 
ergy distribution of the primary (top panel) corresponds 
to that of mid-/late-M dwarfs, with the highest flux in 
the J band, broad water vapor absorption bands cutting 
out the continuum between JH and HK regions and nu¬ 
merous hydride, oxide (FeH, TiO, VO, CO), and alkali 
line (Nai, Ki) absorption features. Also, a Ha emission 
line at 656 nm is detected. The overall appearance of the 
companion spectra (bottom panel) appears to resemble 
an L-type object, but it is significantly redder than field 
L dwarfs. Unlike a typical L dwarf, its flux peaks at 
the H band, which has a triangular continuum shape, 
usually interpreted as a hallmark of low surface gravity 
and youth. In the following, we discuss in detail the 
spectral characteristics, determine the spectral types of 
the binary components, and analyze the gravity-sensitive 
features in the spectra. 


4.1.1. The Primary 

To determine the spectral types, we used our low- 
resolution optical and near-IR spectra. We classified 
the objects in a qualitative manner through visual com¬ 
parison of the spectra with a set of field dwarf spectral 
templates, in the optical and near-IR separately. Subse- 
quently, we u sed the spectral indices established by |Allers| 
& Liu (20131 to designate spectral types in a quantitative 
way, and estimate the gravity class. The final spectral 
types were then assigned by averaging the information 
from both approaches with their corresponding uncer¬ 
tainties. 

In Figure we show the optical (left plot) and 
near-IR (right plot) spectra of the primary overplotted 
with the templates. In the optical, we compared our 
GTC/OSIRIS spectra normalized at 900 nm with a set 
of M dwarf spec tral templates: Cl 406 (M6, [Reid et "ST 


20081), vB 8 (M7, Kirk patrick et al. 19971 and vB 10 (lVi8, 


Henry & Kirkpatrick 


1990|. In this wavelength range 
Hi 


we hnd the best match with an M7-type dwarf, while a 
noticeably poorer fit to templates with one subtype dif¬ 
ference. Thus, in the optical we assign a spectral type 
of M7±0.5. In an analogous manner, we compared the 
near-IR NTT/SofI spectra normalized at 1.6 /rm with M- 
dwarf spectral templates available in the I RTF Spectral 
Librarjrl maintained by Michael Cushing (Rayner et al. 


2009). tJn the right-hand side panel of Figured we in¬ 


clude a comparison of the primary spectra with vB 8 
(M7, 

(M8, 


Cushing et al. 120051 

Rayner et al. 12009 

), LP 412-31 

Rayner et af. [2009^ 

aff smoothed to match resolu- 


tion ot 600 and the M9 standard LHS 2065 (Kirk¬ 


patrick et al. 1991) which we observed using SofI with the 


same conhguration as for VHS 1256-1257. We found the 


® http://irtfweb.ifa.hawaii.edu/-spex/IRTF_Spectral. 
Library/ 


TABLE 3 

Pseudo-equivalent widths (pEW) of lines and doublets 
MEASURED IN THE VLT/UVES, NTT/SOfI AND NOT/ALFOSC 
spectra of THE PRIMARY VHS 1256-1257. 


Line 

Date (UT) 

A (nm) 

pEW (A) 

H, 

2014 Jul 01 

383.5 

-6.06 ± 0.91 

He 

2014 Jul 01 

388.9 

-7.94+ 1.23 

Can K 

2014 Jul 01 

393.4 

-35.63 + 4.42 

Can H-(-H, 

2014 Jul 01 

396.8 + 397.0 

-28.44 + 3.79 

Hs 

2014 Jul 01 

410.2 

-9.49+ 1.78 


2014 Jul 01 

434.1 

-8.43+ 1.87 

H/3 

2014 Jul 01 

486.1 

-7.49+ 0.31 

He 

2014 Apr 22 

656.3 

-3.8+ 0.5 

He 

2014 May 04 

656.3 

-4.0 + 0.5 

He 

2014 Jul 01 

656.3 

-4.04+ 0.05 

He 

2014 Jul 09 

656.3 

-3.74+ 0.19 

He 

2014 Jul 10 

656.3 

-2.80 + 0.09 

Li I 

2014 Jul 01 

670.8 

< 0.03 

K I doublet 

2014 Jul 01 

766.5/769.9 

9.35+0.17/5.94+0.16 

Na I doublet 

2014 Jul 01 

818.3/819.8 

1.80+0.06/2.88+0.06 

Nai 

2014 Mar 12 

1139 

8.2 + 0.8 

K l doublet 

2014 Mar 12 

1169/1177 

3.0+0.6/4.2+0.7 

K l doublet 

2014 Mar 12 

1244 / 1253 

4.0+0.7/3.4+0.6 


best consistency with the M8 standard. Both M7- and 
M9-type standards provided significantly worse matches, 
especially in the J band. We adopt a near-IR spectral 
type of M8 ±0.5, one subtype later than the one desig¬ 
nated from the optical. 

As mentioned above, in parallel to visual comparison of 
the overall spectral morphology with the standards, we 
calculated the spectral indices defined by |Allers & Liu| 
([M^. We used the H2O, H2OD, H2O-I, and H2U-2 
indexes, optimized to avoid spectral features dependent 
on the surface gravity of the object and to have a well- 
established correlation with optical spectral types. The 
H2O index, which measures the slope of the steam ab¬ 
sorption at the blue end of the H band (~1.50-1.57 p.m), 
yields a spectral type of M7.6 ± 0.4. From H2O-I and 
H2O-2 indexes we obtained spectral types of M6.9 ± 1.1 
and M7.9 ± 0.5, respectively. Both from the compari¬ 
son with templates and from the calculated indices, we 
found the typing to be consistent within the uncertain¬ 
ties. As the final spectral type we adopt an M7.5 ± 0.5, 
obtained as the mean of all the types inferred from visual 
comparison and from spectral indexes, weighted by their 
corresponding uncertainties. 

To find out if the primary has characteristics of low 
surface gravity, we inspect the gravity-sensitive spectral 
features like the Nai and Ki lines. In Fig. f^we present 
a comparison of Ki and Nai doublets in the optical at 
around 770 and 819 nm, respectively (top panels of the 
figure), and in the near-IR in a region of the J band (bot¬ 
tom panel). In the optical, we plot in black the UVES 
spectrum of the primary VHS 1256-1257, in red the stan¬ 
dard M7 dwarf vB 8 (using UVES spectra retrieved from 
the ESC Archive) and in blue the M 8 member of the 
Pleiades, Teide 1 (Rebolo et al. 1996). In the near-IR, 
our NTT/SoH spectrum of VHS 1256-1257 (black line) 
is overplotted with the field dwarfs vB 8 (M7), LEHPM 
2-436 (M8) observed by us with the same instrumental 
confi guration, and w ith an M8 at 5-10 Myr, 2M 1207 A 
from Bonnefoy et al. (2014). The spectra were smoothed 
to match the resolution of VHS 1256-1256 data and, in 
the optical, convolved with a Gaussian function to adjust 
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Fig. 6.— Low-resolution GTC/OSIRIS optical (left) and NTT/SofI near-infrared (right) spectra of the companion VHS 1256-1257 b 
plotted with black line, compared with high-gravity field dwarf templates. The sources and references of used templates are described in 
Section 4.1.2. Spectra were normalized at 0.9/rm at optical and at 1.6/rm at near-IR and offset by a constant for display. 



Wavelength (pm) 

Fig. 7.— Sofl/NTT low-resolution near-infrared spectra (0.9-2.4 ^m, ij~600) of the companion VHS 1256-1257 b (plotted in black) 
compared with an L4 it 1 type dwarf SIMP2154-1055 classified as intermediate surface gravity object (INT-G) and with a set of young, 
very low-gravity (VL-G) L dwarfs (red lines). Spectra were normalized at 1.6 fim and shifted by a constant for display. The spectrum of 
VHS 1256-1257 b is best matched with spectrum of the free-floating planet PSO J318.5-22 classified as L7 it 1 type. 


the FWHMs of the lines. 

We find these absorption features of the primary to 
be slightly less intense than in the field counterparts, 
for example, the Nai lines at 818.3 and 819.8 nm have 
pEW of 1.80 ± 0.06 and 2.88 ± 0.06 A, respectively. In 
comparison, from the vB 8 spectra we measure pEW of 
1.98 ± 0.07 and 3.40 ± 0.06 A. The same lines in the 120 
Myr old M8-type counterpart Teide 1 are in turn less in¬ 
tense with pEW of 1.66±0.09 and 2.66±0.10 A. Because 


the Nai subordinate lines are well known to be highly 
sensitive to atmospheric pressure at cool temperatures, 
this indicates that the primary likely has an intermedi¬ 
ate surface gravity that lies between that of Pleiades and 
field M8-type dwarfs. We note, however, that all these 
comparisons are based on the assigned spectral types and 
assume that both targets and reference sources have the 
same metallicity. Eurthermore, the intensity of the al¬ 
kali lines also depend on the effective temperature (or 
spectral type). Given a half-subtype uncertainty in our 
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typing of VHS 1256-1257 A, we found that the surface 
gravity of this object may be consistent with field ob¬ 
jects considering an earlier type. 

In Table we provide the pEW of Nai, Ki, and other 
lines measured in the optical and near-IR spectra of 
the primary. In order to quantify the gravity class we 
have calculated four spectral indices: FeH^, YOz, Kij, 
and i7-continuum, identified to be related to gravity- 
dependent features and optimiz ed to be employed w ith 
low-resolution near-IR spectra (Allers & Liu 2013 and 
references therein). Based on these indicators we consid¬ 
ered field gravity (FLD-G), intermediate gravity (INT- 
G), or very low gravity ( VL-G) of an object f ollowing the 
classification scheme of Allers & Liu (2013). For FeHz, 
which measures the depth ot the Fell feature at 0.99 /xm, 
we obtained a value of 1.13 ± 0.08. For the VO^ index 
we assign a score of “n” since it is not gravity sensitive 
in this range of spectral types. For the Kij index, which 
measures the depth of potassium absorption at 1.244 and 
1.253 ^m, we got a value of 1.08 ± 0.(17 and for the H- 
cont index, which measures the slope of the H band we 
got 0.98 ± 0.24. For an M7.5 spectral type these index 
values correspond to a surface gravity intermediate be¬ 
tween that of counterparts classified as very low gravity 
and of the old, held population. 

4.1.2. The Secondary 

Our GTC/OSIRIS optical spectrum of the compan¬ 
ion (left plot of Fig. is of modest signal-to-noise 
ratio, and apart from the general spectral energy dis¬ 
tribution only few spectral features are noticeable, like 
the strong Ki resonance doublet and water vapor at 
0.92 /xm. The lack of oxide features (TiO, VO) at 
the resolution of the visible spectrum suggest a type 
later than mid-L. On the left panel of Fig. |^we show 
the spectra of the companion (black line) normalized at 
0.9 /xm and overplotted with obj ects of spectral t ypes 
L7 (2MASS J21522609-H0937 575, iReid et al.|[2 00^, L8 


(2MASS J03400942-6724051, Reid et al. '2UD8) and TO 


(2MASS J04234858-0414035,|Gruz et al.|2003|. The tem¬ 
plate spectra were retrieved from the Ultracool RIZzo 
Spectral Library. Gomparison with templates of held 
dwarfs yields a spectral type of L8±2. 

From the comparison of the NTT/Sofl near-IR spec¬ 
trum of the companion VHS 1256-1257 b to L-type ob¬ 
jects with well-determined spectral types, we hnd that 
the overall 1.0-2.4 /xm spectral energy distribution does 
not match any early to late high-gravity held L dwarf. In 
the right-hand side plot of Fig. we show a comparison 
of the companion near-IR spectra normalized at 1.6 /xm 
(black line) with an L3-type o bject DENIS-P J1058.7- 


1548 (Bu rgasser et al. 
005902. L) dCushing et ah 


5319A (Luhman 2013). 


2010), L5 SDSS J053951.99- 


mS) and L8 ± 1 WISE 1049- 
Ihe spectra of these L dwarf 
templates were obtained by us using SofI instrument with 
the same conhguration as for VHS 1256-1257 and were 
reduced in the same way. In the J band, VHS 1256- 
1257 b is relatively less luminous than normal held L 
dwarfs and contrarily is signihcantly brighter in the K 
band. Its iL-band continuum unlike the held L dwarfs 
has a distinctive triangular shape, interpreted as a sig¬ 
nature of low surface gravity and youth. To measure 
the near-IR spectral type, we tried to hnd best match¬ 
ing held dwarf standards from the IRTF library, at the 


J (1.07-1.40 /xm) and K (1.90- 2.20 /xm) windows s epa- 
rately, following the approach of Allers & Liu (2013). In 
this wavelength regions the shape of spectra continuum 
is expected to have a lower dependence on gravity. We 
estimate, in the J and K bands, respectively, a spectral 
type of L8 ± 1 and L5 ± 2. In the K band we assign an 
uncertainty of two subtypes because several standards 
provide similarly good hts. 

Along with visual comparison, we calculated the H2OD 
index, well dehned up to the range of late-L dwarfs. 
The other indexes are not valid for spectral types later 
than L5. The H2OD index value give a spectral type 
of L5.8 ± 0.8. In analogous manner as for the primary, 
from a mean of the types inferred through visual compar¬ 
ison and from the H2OD spectral index, weighted by the 
corresponding errors, we obtained an L7.0 type, which 
we adopt here as the spectral type of the companion, 
with an uncertainty of 1.5 subtype. Among the gravity- 
sensitive spectral indices, only the i7-continuum index is 
well defined up to L7 type objects. For VHS 1256-1257 b 
iJ-cont i ndicates a very low g ravity (VL-G). However, as 
noted by Allers & Liu (2013), the shape of the i7-band 
continuum is not the most reliable indicator, and can 
provide rather only a hint of low gravity. It should be 
used in combination with other gravity-sensitive indices, 
but for the spectral type in question none of them have 
yet been defined. The Ki lines at 1.17 /xm (pEW < 5 A) 
and at 1.25 /xm (pEW < 3 A) appear weaker in the spec¬ 
trum of the secondary than in the field mid- and late-L 
dwarfs (see Fig. |^. For example, in the spectrum of the 
L8 dwarf WTS'F 1049-5319 A obtained with the same in¬ 
strument and setup we measured a pEW of the K i lines 
of 5.7±0.8 and 6.8±0.7 A at 1.169 and 1.177 /xm and 
3.4±0.4 and 6.1±0.4 A at 1.244 and 1.253 /xm, respec¬ 
tively. This also indicates that the secondary has a lower 
gravity compared to field counterparts (?). 

In addition to the comparison with held L-type tem¬ 
plates, we have also compared the near-IR spectra of 
VHS 1256-1257 b with a set of known young low-gravity 
mid- and late-L dwarfs. In Figure we overplot the 
companion spectra with SIMP J2154-1055, class ihed as 
an L4 ±1/3 with intermediate surface gravity (Gagne 
let al.||2014|), a ~5 Mj^p pl anetary mass companion 


-3932 b at 8 Myr (Ghauvin et al. 


2005 


Pa- 


assi- 


tience et al. 2010), 2MASS J22443 167±21MT433~c 
hed as very low-gravity L6± 1 type (McLean et al.|2003l 


and with the free-hoating planetary mass object PSO 
J318.5-22 (Liu et al. 2013). In general, the spectrum of 
the companion hts rnuch better to the spectra of young, 
very low-gravity L dwarfs than to the spectra of held ob¬ 
jects. We hnd the best match with PSO J318.5-22 which 
is classihed as an L7 ± 1, thus providing a further sup¬ 
port of our spectral type classihcation. PSO J318.5-22 is 
a 6.5 Mjup object, possible member of the 12 Myr mov¬ 
ing group /3 Pictoris. Because of this strong resemblance, 
we state that the companion can be assigned to the very 
low-gravity (VL-G) class. 

4.2. Distance Estimates 
4.2.1. Spectrophotometric Distance 

Having classihed the spectral type of the primary, we 
used the JHKg photometry to estimate the distance of 
the system. We used the 2MASS catalog photometry. 
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since at these magnitudes the VHS measurements start 
to get beyond the linear regime of the detector and could 
be uncertain. We employed the mean absolute magni¬ 
tudes as a fu nction of spectral type determined by [Dupuis 


&: Liu|(|2012|) from a set including 8 M7.0, 9 M7.5, and 11 
M 8 .U type dwarfs with precise parallactic distance deter¬ 
mination. Considering an M7.5±0.5 type of the primary, 
and assuming that it is a single object, we inferred spec¬ 
troscopic distances of 12 . 0 ^ 3 ®, 12A^fl and 12 . 2 ^ 3 ;^ pc, 
using J, H, and Kg photometry, respectively. Taking 
the average of three bands, we estimate a distance of 
12 . 2 l 3'5 pc. The quoted errors account for uncertainties 
of the photometry, spectral type, and the intrinsic dis¬ 
persion in absolute magnitudes for field dwarfs at a given 
spectral type. 

In this case however, the estimation that employs the 
M dwarfs found in the field should be considered with 
caution, since the field objects have ages typically of 
more than 1 Gyr, and for VHS 1256-1257 we expect an 
age below 300 Myr (see Section 4.5). From the studies 
of young star clusters, it is known that the M dwarfs 


respect to their field counterparts (Bihain et al. 

2006 

Jameson et al. 

2008 

/apatero Usorio et al. 

2014c 

1 . For 


that reason we obtained an alternative estimation using 


Pleiades cluster, compiled by Bihain et al. (|2010b The 

Pleiades has an age of 120 ± 10 Myr (Basri et al. 11996 

Martin et al.|1998 

Stauffer et al.|1998D and is located at 

a distance of 133.5 

pc (|Soderblom et al.|2005 

Melis et al. 

2014 

1. Averaging the estimates obtained in 

J, Jd, and 


Kg yields a distance of 14.9_3;2 pc for VHS 1256-1257 if 
this binary had the age of the Pleiades. 


4.2.2. Trigonometric Parallax 

To derive the trigonometric parallax and proper mo¬ 
tion of the pair we used an area of 5' x 5' around the 
target and the images listed in Table except for the 
2MASS. We took advantage of the avSability of these 
images for a first determination of a trigonometric dis¬ 
tance to the system. In Table|^we provide the observing 
epoch (including the Julian Date), telescope, and imag¬ 
ing instrument, filter, the average pixel projection onto 
the sky, and the mean seeing of the images. All data were 
properly reduced following standard steps for the optical 
and near-IR wavelengths as explained in Section 3. We 
selected the VISTA/VIRGAM V-band observations of 
2011 July as the fundamental frame to which all other im¬ 
ages are compared. Using the daofind command within 
IRAF we identified all sources with photon peaks with 
detection above 6 cr, where a stands for the noise of the 
background, and FWHM resembles that of unresolved 
objects (i.e., extended sources were mostly avoided). In 
addition, we ensured that the detected sources lied within 
the linear regime of the detectors response (with the only 
exception of the primary component of the binary). The 
number of sources identified per image in common with 
the reference frame ranged from 25 to ^40. The centroids 
of detected objects were computed by estimating the x 
and y pixel positions of the best htting one-dimensional 
Gaussian functions in each axis; typical associated errors 
are about 3%-5% of a pixel or better. 

Pixel coordinates were transformed between different 



Fig. 8.— Relative apparent trajectory of the components of 
VHS 1256-1257 in right ascension (a) and declination (5) coordi¬ 
nates for the time interval 2011 July (reference epoch) through 2014 
December. Each component is plotted separately as indicated in 
the figure legend, and the averaged values per observing month are 
shown with the red solid circles. The solid line indicates the best 
model fit. North is up and east is to the left. 


epochs using the GEOMAP routine within iRAF, which ap¬ 
plied a polynomial of the third and fourth order in x and 
y and computed linear terms and distortions terms sep¬ 
arately. The linear term included an x and y shift and 
an X and y scale factor, a rotation, and a skew. The dis¬ 
tortion surface term consisted of a polynomial fit to the 
residuals of the linear term. The {x, y) astrometric trans¬ 
formation between observing epochs and the reference 
epoch was an iterative step, which included the rejection 
of objects deviating by more than 1.5-2 cr, where a cor¬ 
responds to the dispersion of the transformation. Typ¬ 
ical coordinates transformation dispersions ranged from 
±0.05 to ±0.09 pixels; this is typically ±15-25 mas. The 
centroids of the two components were calculated with an 
uncertainty similar to that of the astrometric reference 
sources. Errors associated with the astrometry are domi¬ 
nated by the errors in the coordinate transformations for 
each axis (R.A., deck), which include the dispersion of 
all reference sources. We therefore assigned the disper¬ 
sion of the coordinates transformation solutions to the 
errors associated with the differential astrometry of our 
targets listed in Table 4. We did not apply any correction 
for differential chromatic refraction since all of our data 
were taken at red wavelengths, where the corrections are 
smaller than the quoted astrometric uncertainty. The 
relative {dx, dy) astrometry (conveniently derotated) was 
converted into da cos 6 and dJ using the corresponding 
plate scales. The apparent trajectory of each component, 
which depends on proper motion (/Tq., Ts) and parallax 
(tt), is shown in Figure]^ and can be modeled with the 
following equations: 

da =/i„ (t-to) ±7r (/“-/“) (1) 

d5 = ys {t-to)+T: {f! - ft), ( 2 ) 

where t stands for time, the subscript o indicates the 
reference epoch, and /“ and refer to the parallax fac- 
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tors in right ascension (a) and declination (J), respec¬ 
tively. In our study, all the astrometric quantities are 
given in mas and the times t and to are measured in 
Julian Days. The parallax fac tors we r e com puted by fol¬ 
lowing the equations given in Green] ( 1985] ) and obtain¬ 
ing the Earth barycenter from the DE405 Ephemeris]^ 
We applied the least-squares fitting method to the set 
of Equations and to derive the parallax and the 
proper motion of the system. The two components 
were fit simultaneously. The best-fit solution yielded 
/Tq, cosS = —277± 5 masyr“^, fis = —203± 12.5 masyr“^ 
giving a total proper motion of ^ = 344±13 masyr“^, po¬ 
sition angle of the proper motion vector of 233.7° ±2.5°, 
and a relative parallax 7 r = 77.79 ±6.4 mas. It is de¬ 
picted with a solid curve in Figure The amplitude 
and position angle of the proper motion we deter mined 
is in agreement with the measurement given by Dea- 


con & Hamblyj (2007). From the fitting performed for 
each pair member individually we obtained proper mo¬ 
tion of fia cos <5 = —281.5 ± 5.3 masyr“^, fig = —205.5 ± 
15.2 masyr“^ for the primary and /Tq: cos <5 = —275.4 ± 
5.3 masyr“^, ^5 = —198.4 ± 15.2 masyr“^ for the com¬ 
panion. The above values are listed in Table We 
have also tested the WISE, 2MASS, and available earlier 
epoch data (DENIS, USNO-Bl.O, DSS) in the proper mo¬ 
tion and parallax determination, but despite the longer 
baseline they introduced larger errors to the fitted values. 
Hence, for the determination of the parallax and proper 
motion we excluded these data, given that their associ¬ 
ated uncertainties are at least one order of magnitude 
higher than those of the recent data. We remark that 
given the time coverage of our data, which span 3.0447 
yr, the proper motion is more precisely determined than 
the parallax. Additional images, in particular taken on 
months that were not covered yet, are required for de¬ 
creasing the uncertainty in the parallax measurement. 

We did not correct our proper motion for the motions 
of the stars used as a reference because we assumed that 
these (small) motions are randomly orientated with a 
negligible net effect. However, we did apply the correc¬ 
tion for converting the relative parallax into the abso¬ 
lute parallax. This takes into account that the reference 
stars are located at finite distances, which diminish part 
of the true parallax of our targets. We followed a proce¬ 
dur e similar to the one described in Faherty et al. (2012 1 
and jZ apatero Osorio et al.j (|2014a|) Using the ziVlAyS 
colors ( jSkrutskie et al.|2006 ) of the reference objects, we 
obtained their photonietric distances by assuming that 
all of them are main sequence stars. We adopted the 
co lor-bolometric c orrection-spectral type rel ations given 


Johnson (1966) for BAFGK stars and in {Kirkpatrick 
et al. p99;j| tor late-K and M stars. The defined rela- 
tions are valid for colors in the interval J — Kg = —0.2 to 
1.53 mag. We adopted the mode of the distribution of 
reference objects distances as the correction to be added 
to the relative parallax that comes directly from our fit 
to obtain the absolute parallax. The absolute parallax is 
7 r = 78.79 ±6.4 mas, which translates into a distance of 
12.7 ± 1.0 pc, consistent with the value from the spec- 
trophotometric estimates using field counterparts. This 
suggests that the primary component is not an equal- 
mass binary. 


http://ssd.jpl.nasa.gov 


4.3. Angular Separations and Orbital Motion 

On the available images where both components were 
detected and the primary did not saturate we measured 
the projected angular separation (p) and position angle 
(d) of the companion using the centroid positions of the 
two sources in each of the images transformed into the 
VHS a, 6 coordinates as explained above. The deter¬ 
mined values are given in Table The angular separa¬ 
tions and position angles of VHS1256-1257 in the most 
recent epoch images over the last three years are con¬ 
sistent within 6p < 30 mas and 69 < 0.04 deg, respec¬ 
tively. Both components share the same proper motion, 
and given the low probability to find a very red young L 
dwarf (only a few tens found in the whole sky area) at 
this very short angular separation of 8706, we conclude 
that both objects are gravitationally related. These re¬ 
sults also indicate that we do not detect an orbital motion 
higher than 10 masyr“^ (at the Icr level). 

We expect an orbital period of about 3900 yr from the 
estimated masses of the primary and companion (Sec¬ 
tion 4.7) and a physical projected separation of ^102 AU. 
Assuming a circular orbit, this implies that the displace¬ 
ment caused by orbital motion would be of 4-13 mas yr“^ 
depending on the orientation of the orbit, which is con¬ 
sistent with our results. The angular separation of the 
pair in the 2MASS images is lower, although consistent 
within 3cr. These differences ((y.'6±072) can not be ex¬ 
plained by the orbital motion of the companion, since we 
expect no more than ~ 13 masyr“^ for a face-on circu¬ 
lar orbit at the physical separation found, and may be 
attributed to systematic errors in the determination of 
the centroids in the 2MASS data probably due to small 
separation of the components and low resolution of the 
images. 


4.4. Kinematics 

To measure the heliocentric radial velocity of 
VHS 1256-1257 we employed the high spectral res¬ 
olution UVES data (mean Modified Julian Date, 
MJD = 56838.3530) and the cross-correlation method 
against the M 6 V star GJ 406, which has a known , con¬ 
stant radial vel ocity of = 19.5±0.1 kms“^ (Nide- 
ver et al. 2002 ) and a small projected rotational ve - 
locity of vs'mi < 3 kms“^ (Reiners & Basri 2010). 
GJ406 was also observed with the VLT/UVES instru¬ 
ment (spectral resolution of ~45,000) on 2009 March 14 
(MJD = 54904.0839). We downloaded the reduced UVES 
spectrum of GJ 406 from the ESO data archive (program 
082.D-0953). Although GJ 406 was observed with a non¬ 
standard i nstrumental setup in ter ms of wavelength cov¬ 


erage (see Fuhrmeister et al. 


the M 6 V-type star and VHS 


2011 ), the overlap between 


256-1257 data is significant 


and secures a reliable velocity determination. The tel¬ 
luric contribution was removed from the data of GJ 406 
in the same manner as we did for VHS 1256-1257 (see 
Section 3.4). 

The telluric-free UVES spectra of VHS 1256-1257 and 
GJ 406 were cross-correlated using the task FXCOR within 
IRAF. We fit a Gaussian function to the peak of the 
cross-correlation distribution. The resulting relative dis¬ 
placement was corrected for the lunar, diurnal, and an¬ 
nual velocities to obtain the heliocentric radial velocity 
of VHS 1256-1257. All UVES spectra were calibrated in 
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wavelength using ThAr arc lines taken some time after 
the science observations (Section 3.4). The typical sta¬ 
bility of UVES along the spectral axis is on the order of 
a pixel over the course of a night; therefore, deviations 
of a kms“^ can be expected. To correct our velocity 
determination for this effect, we employed the absorp¬ 
tion lines from the Earth’s atmosphere present in the 
data, which provide a reliable refere nce frame within ±90 
ms“^ (e.g., Griffin & Griffin||l973 ). We used the UVES 
spectra of VHS 1256-1257 and GJ406 before subtrac¬ 
tion of the sharp telluric lines to determine the instru¬ 
mental shifts of the wavelength solution for each object. 
We found them to be ±0.62 ±0.20 kms“^ (VHS 1256- 
1257) and -0.14±0.20 kms"! (VHS 1256-1257 versus 
GJ406). The final corrected heliocentric radial velocity 
of VHS 1256-1257 is 1 ;^ = —1.4±5.0 kms“^, where the 
error bar accounts for the uncertainties due to the cross¬ 
correlation procedure, the corrections of the wavelength 
reference frames, and the error associated with the ve¬ 
locity of the M6V standard star. 

We checked this heliocentric velocity using the line cen¬ 
troid method, which is independent of any reference star. 
By measuring the centroids of alkali and other metal¬ 
lic lines in the three UVES spectra of VHS 1256-1257 
and considering the correction for the wavelength refer¬ 
ence frame and lunar, diurnal, and annual velocities as 
previously discussed, we found w/j = ±4.1±4.0 kms“^. 
Albeit having a much lower spectral resolution, the 
NOT/ALFOSC spectra of VHS 1256-1257 also yielded 
radial velocities that agree with the UVES result at 
the ItT level (where a stands for the velocity uncer¬ 
tainty): = ±0.1 ± 7.0 kms“^ (MJD = 56769.0443) and 

±14.1 ± 10.0 kms-i (MJD = 56781.9229). The ALFOSG 
velocities were derived using the cross-correlation tech- 
nique; the M4.5V star GJ388 {vh = ±12.6215 kms~^ 


Bonfils et al. (2013); Vh = ±12.453 kms Chubak et al. 


2012p acted as the velocity reference source because its 


spectra were acquired with exactly the same instrumen¬ 
tal configuration and observing dates as our target (Sec¬ 
tion 3.3). G J 388 was found to b e a radial velocity vari¬ 
able star by Bonfils et al. (2013). However, the ampli¬ 
tude of the variations is 0.25 kms“^, which is about >30 
times smaller than the velocity precision we were able 
to achieve with the ALFOSG data. No obvious velocity 
change larger than ^10 kms“^ is observed in VHS 1256- 
1257 over the 69 day interval of spectroscopic observa¬ 
tions. 

Having the measurements of proper motion, parallax, 
and radial velocity we can determine the three compo¬ 
nents of the Galactic space veloc ity, U, V, and W. We 
applie d the formulas presented by Johnson fc Soderblom 
(]1987P to derive the UVW velocities listed in Table [2 
and displayed in Figure]^ The uncertainties associatea 
to all three Galactic velocities come from the proper mo¬ 
tion, parallax, and radial velocity error bars. Figure 
also illustrates the ellipsoids corresponding to well-known 
young stellar movin g groups of the solar neigh bor hood 


(data compiled from Zuckerman & Song 2004 and Tor- 
res et al.]|2008 |. 


As seen m Figure the space velocities of VHS 1256- 
1257 fits with the /3Pic Group and is consistent within 
error bars with the galactic velocities of Taurus-Auriga, 
TWHya, and the Local Association groups. The stars 






Fig. 9.— Galactic space velocities of VHS 1256-1257 (red dots) 
with overplotted ellipsoids of known young star associations and 
moving groups. Errors incorporate uncertainties in the proper mo¬ 
tion, parallactic distance, and radial velocity. Galactocentric U 
velocity is positive toward the Galactic center. 


in P Pic, Taurus-Auriga, and TW Hya are younger than 
20 Myr, and this age is incompatible with the non¬ 
detection of lithium in the atmosphere of the primary 
(see next section). Hence the kinematic and spectro¬ 
scopic properties of VHS 1256-1257 is only compatible 
with the membership to the Local Association Group, a 
coher ent kinemati c stream of young stars (all below 300 
Myr, Eggen 1992) with constituent clusters and associ- 
ations such as the Pleiades (120 Myr), a Persei (50-80 
Myr), and IC 2602 (~70 Myr). Membership in stellar 
moving groups is commonly used to constrain or confirm 
the ages especially of young (< 1 Gyr) objects. We dis¬ 
cuss in more detail the age of VHS 1256-1257 system in 
the following section. 


4.5. Activity, Lithium and Age 

T he high -resolution optical spectra of the primary 
(Fig.[W 11) reveals emission features recognized as indi¬ 
cators of chromospheric activity in very low-mass stars. 
We detect the Balmer series emission lines of hydrogen, 
from H(j up to H^, and single ionized calcium Ga ii H and 
K emission lines at 396.8 and 393.4 nm. Measurements 
of the lines pEW are given in Table The pEW of H^ 
does not show a strong variability in timescales of days 
and months, but show small variations of the order of 


11 


1-2 A along the observations, as can be seen in Fig. 
and Table These activity indicators of the primary 
are consistent with dwarfs of similar spectral types in 
the field and in intermediate age clusters l ike the Hyades 
(Reid & Hawley 1999 West et al. 2011), and are less 
active than those of younger open clusters like Pleiades 
( Stauffer et al.||1998 |. 

In the VL'iy (j VES spectrum of VHS 1256-1257, we did 
not detect the Lii resonance doublet at 670.8 nm, impos¬ 
ing an upper limit of 30 mA in the pEW of this atomic 
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Fig. 10.— VLT/UVES spectra of the primary VHS 1256-1257 covering 373—495 and 565—950 nm wavelength range at a resolution of 
/?~40 000 (25 mApix”^). The blue part of the spectrum (left panel) was smoothed for the display by a factor of 33, degrading the 
resolution to 10 000. The spectrum has been corrected for telluric absorptions using the ESO Molecfit software. Certain features 

indicating strong magnetic activity in the chromosphere are visible: the Balmer series emission lines from up to Hrj, singly ionized 
calcium H and K lines at 396.85 and 393.37 nm, respectively. 


line (see the right panel of Fig.[IT|). This is much lower 
than the expected value for a fulTpreservation of this el- 
eme nt in the atmosphere of late-M dwarfs of EW = 0.5 - 
lA (Stauffer et al. 1998 Zapatero Osorio et al. 2002). 


The Li element is rapidly destroyed in the interior of 
stars, on timescales shorter than ~ 150 Myr, and in mas¬ 
sive brown dwarfs on timescales of a few Gyr. Brown 
dwarfs with masses lower than 0.055-0.060 Mq do not 
burn this element in their interiors because their cen- 


tral temperature is not high enough to 

rroduce this fu- 

sion reaction ( 

Ushomirsky et al. 

1998 

Ghabrier et al. 


pose a lower limit in the mass of the primary of 0.055- 
0.060 Mq and also imposes a lower limit to the age of 
the system, since objects of similar spectral type in the 
Pleiades cluster (age^ 120 Myr ) have fully preserved this 
element (St auffer et al. 1998|)' Theoretical evolution¬ 
ary m odels (Ichabrier fc Baratte 1997 Chabrier et al. 
2000b) predict that objects with ettective temperature 
'1 eff 2600 K (which is the corresponding effective tem¬ 
perature of an M7.5) have preserved their Li content in 
timescales of less than 150 Myr, but have destroyed it for 
several orders of magnitude at ages larger than 200 Myr. 

The gravity-dependent spectral features of the primary 
like Na and K alkaline lines are also consistent with inter¬ 
mediate grayities_between_thos£_of_theKemdesand_fiUd 


dwarfs ([Go rlova et al. |2003[ lAllers fc Liu|2013[|Bonnefoy| 
et al .|20Tir Galvez-Urtiz et al.|2014|) Un the other hand, 
second; 




the secondary shows spectral features, which are sign¬ 
posts of low gravity and youth like the sharp triangular 
shape of the H band. Some authors attribute thes e fea¬ 
tures to ob.iects with ages young er than 150 Myr (Gruz 
et al.||2009| Faherty et al.||2013a), which is in contradic- 
tion with the lower limit to the age of VHS 1256-1257, 
given by the non-detection of Li in the primary. Other 
authors sugge st a more conservative range of ages up to 
120-500Myr ( Zapatero Osorio et al.|2014a ), which are in 
good agreement with the age constrain tor this system. 

The proper motion, radial velocity, and parallactic dis¬ 
tance of the primary allow us to determine its galac¬ 
tic kinematic. The galactic velocities UVW of the 
VHS 1256-1257 indicate that the system probably be¬ 
longs to the Local Associat ion, whose membe rs have esti¬ 
mated ages of 10-300 Myr (Asiain et al. 1999). In conclu¬ 



Fig. 11.— Close-up of regions of the VLT/UVES spectrum of the 
primary VHS 1256-1257, showing, in the left panel the Ha emission 
line at 656.3 nm and in the right panel the non-detection of lithium 
absorption at 670.8 nm. Flux uncertainties are plotted in gray. 


sion, based on the absence of lithium in the primary and 
the likely membership to the Local Association, we adopt 
a range of 150-300 M yr for the age of the VHS 1256-1257 
system. 


4.6. xs Versus J — Kg Color-Magnitude Diagrams 

Average near-IR J— Kg color of L7 dwarfs is 1.75±0.26 
(Faherty et al. 2013b). With J — Kg = 2.47±0.03 mag, 
VHS 1256-1257 b is among the reddest known L dwarfs, 
next to 2M1207 b and PSO J318.5-22. These extreme 
colors are thought to originate from enhanced photo- 
spheric dust and broad-band changes in the spectral en¬ 
ergy distribution induced by the low gravity of young 
objects. In Figure we compare the absolute J and 
Kg magnitudes versus J — Kg colors of VHS 1256-1257 
with held mid-late-M, L and T dwarfs with measured 


parallaxes compiled by Dupuy & Liu (2012) and with 
several known substellar objects and giant planet com¬ 
panions. Based on parallactic distance we derived an ab¬ 
solute magnitude of the companion Mj = 16.45 ± 0.30, 
which is about 2.2 mag fainter than the late-L held coun¬ 
terparts. In the Kg band it is roughly one magnitude be¬ 
low the LT sequence. The location of VHS 1256-1257 b 
on the Mj and Mks versus J — Kg color magnitude dia¬ 
grams coincides with the young planetary mass compan¬ 
ion 2M 1207 b and is almost the same as the 5-11 Mjup 
planet HR 8799 b. 
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Fig. 12.— Mj and Mj^s vs. J — Ks color-ma gnitude diagrams com paring the two components of VHS 1256-1257 with field M, L, and T 
dwarfs with known parallax measurements fromi Dupuv &: Liu | (|2012[l. known y oung substellar objects from the compilation of |Bowle r et al.| 
st-massive Pleiades members jzapatero Usorio et al.|2014b[|. All photometry was converted to the MKO filter system. 


l|2013|l and the least-n 


4.7. Luminosity, Mass, and Effective Temperature 

To derive the bolometric luminosities of the primary 
and the companion, we used the near-IR photometry 
and the parallactic distance. For the primary, which 
does not deviate strongly from the held sequence, we 
have obtained the bolometric magnitude applying the 
corresponding BCj and BCic b olometric corrections de¬ 
ter mined for held obje cts, from Golimowski et al. (20041 
and Dahn et al. (2002|). These BGs combined witn abso¬ 
lute magnitudes and Mboi = 4.73 mag for the Sun yield 
the luminosity of log(Lboi/T 0 ) = —3.f4 ± 0.10 dex. The 
error accounts for the uncertainties in distance, photom¬ 
etry, and bolometric correction. Since the near-IR pho¬ 
tometry of the companion differs signihcantly from the 
“normal” held L dwarfs, the bolometric corrections de- 


termined for held ultracool dwarfs are not valid ( Luhman 


2012 

Faherty et al. 

2012 

Zapatero Osorio et al.||2014a) 

We applied bolomel 

ric corrections (JHKs, 2MASS sys 


tern) derived from the measurements of PSO J318.5-22 


(Liu et al. 2013), which shows a strong similarity to 


VHS 1256-1257 b. Taking the mean value from the three 
bands, we obtained log(Lboi/L 0 ) = —5.05 ± 0.22 dex. 

We used the calculated bolometric luminosities to infer 
the range of possible masses, effective temperatures, and 
gravities (log g) from the evolutionary models, adopting 
an age of the system of 150-3 00 Myr. We have used the 


Lyon group models BT-Settl (Baraffe et al. 1998 2003 


Ghabrier et al.|2000b ) with the Gattau et al. (2011) solar 


abundances and th e ‘hot start’ cloudy atmo sphere evo¬ 
lutionary tracks of Saumon & Marley (2008), with sed¬ 
imentation parameter }sed = '^- Both models give con¬ 


sistent values of the derived parameters (masses. Teg, 
log g) within their corresponding uncertainties. For the 
primary, we obtained a mass of 731^° ALjup, close to the 
boundary between stars and brown dwarfs regime which 
starts at masses below ^ 75 Mjup for solar metallici- 
ties. The effective temperature and log g found from the 
models are 2620 ± 140 K and 5.05 ± 0.10 dex, respec¬ 
tively. This Teff is consistent with typical temperatures 
of field M7.5 ± 0.5 dwarfs being in the ^25 00-2600 K 
range ( ILeggett et al. 2000 Golimowski et al.||2004 ). The 
derived gravity is slightly lower but similar to the ex¬ 
pected gravity for field dwarfs of similar spectral types 
(log 5 = 5.0-5.5). This is consistent with what we have 
found spectroscopically. 

For the companion, we obtained a mass of 
11 . 2 l® g Mjup, indicating that the object is near the 
mass limit at which the onset of deuterium fusion oc¬ 
curs. Given the uncertainty in the mass determination, it 
is currently unclear whether the object is above or below 
it. The effective temperature and log g of the compan¬ 
ion obtained from the evolutionary models are 880^^^° 
K and 4.241 q'^q dex. A typical effective temperature of 


field L 6 -L 8 dwarf is between 1600 and 1300 K (Leggett 


et al. 2000 Golimowski et al. 2004 Vrba et al. 2003 ). 'Lhe 
effective temperature ol the companion indicated by the 
models is about 500 K lower than that and corresponds 
rather to a field mid-T dwarf. A similar discrepancy 
between temperatures predicted by evolutionary models 
and those estimated from spectral classification has been 
fou nd in other you ng substellar objects: PSO J318.5- 
22 (iLiu et al.||2C)13|), 2MASS 0122-2439 B (iBowler et al. 
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0.02 Ms„„ 


Fig. 13. — Luminosity and aee of the two componen ts of VHS 1256-1257 compared with evolutionary tracks from th e solar abundance. 
cloudy atm osphere models of [Saumon &: Marlevl l|2008[l with /ged = 2 and from the BT-Settl models of the Lyon group (|Allard & Barman] 

H . The |Saumon &; Marleyl (|2UU8| models are plotted in color, with 0.002-0.010 Mq tracks in green (spaced by 0.002 Mq), O.Oli- 
Mq tracks in yellow witti a U.UUl Mq step, and 0.02-0.08 Mq in blue, with 0.01 Mq increments. The light gray curves correspond 
to BT-Settl models . For comparison purposes we overplotted the population of known young substellar objects , from the compilat ions by 
Bowler et a l.ll|2013 1. iLiu et al.l|20]^ and updated with recently discovered objects GLf Psc b ([Nai^et al.|2014f and HD 106906 b ||Bailey| 
et al. 1201411. lire Bi-Betti isomass lines yield slightly higher luminosities at masses > O.OSOMq. 


2013), HN Pe{ 

B (Luhman et al. 

to 

o 

o 

), 2MASS 1207 b 

1 Skemer et al. 

2011 Barman et a 

Ipor 

lb). Moreover, at 


pect to see the methane absorption bands in the near-IR 
spectra of the companion. The absence of methane was 
observed in other similar objects, like for example the HR 
8799 planets. On the contrary, if the actual temperature 
is higher, like that corresponding to “normal” field L 6 -L 8 
dwarfs, this would lead to unrealistically small radius of 
the object (~0.5 Rjup), given its bolometric luminosity. 

We compare in Figure [T^ the luminosity of the com¬ 
panion at the ad opted range of ages with the ISaumon 


&: Marley (2008) and BT-Settl models (Allard & Bar- 


man 


2004 1 and with presently known yo ung substellar 


objects, compiled by Bowler et al. (2013) updated with 
recently reported y oung 'I'-dwart companion GU PSc b 
(Naud et al. 2014) and a pla netary mass companion 
HU 106906 b (Bailey et ah' 2014). Interpretation of evolu¬ 
tionary tracks m this specific region has to be considered 
with caution. At that range of ages and luminosities the 
onset of deuterium burning causes overlapping of isomass 
tracks and less massive object can be slightly more lu¬ 
minous than a more massive one. This may introduce 
ambiguities in the determination of masses and produce 
a relatively large uncertainty in the upper limit of the 
derived mass of the companion. This object for its loca¬ 
tion in the HR diagram and its close distanc e is an ideal 


case for the application of th e deuterium test (Bejar et al. 


1999 Chabrier et al.|[2000a ). 


5. SUMMARY AND FINAL REMARKS 

We have identified an unusually red (J — Ars = 2.47 
mag) common proper motion L7±1.5 type companion lo¬ 
cated at 8"06 ± (/.'OS (~ 102 AU) of an M7.5 ± 0.5 dwarf. 
The near-IR spectrum of the secondary shows a peaked 
triangular shape of the iJ-band continuum, and other 
spectral features recognized as hallmarks of low surface 
gravity and youth. The optical spectrum of the primary 
shows no Li I at 670.82 nm at a limit of pEW < 30 mA. 
We have determined a parallactic distance of the sys¬ 
tem of 12.7 ± 1.0 pc, which is in agreement with the 
distance estimated from spectral type and photometry 
of the primary. From the proper motion, distance, and 
radial velocity of the VHS 1256-1257 we obtained the 
galactic velocities of the primary, which indicate that 
the system likely belongs to the Local Association. The 
non-detection of lithium and the kinematics of the pri¬ 
mary allowed us to constrain the age of the system in 
the range of 150-300 Myr. From near-IR photometry 
and bolometric corrections we estimate luminosities of 
log(Lboi/A 0 ) of -3.14 ± 0.10 and -5.05 ± 0.22 dex for 
the primary and secondary, respectively. By comparison 
with theoretical evolutionary models we derived a mass 
of 73^11 Afjup for the primary, at around the substellar 
mass limit and 11.2'^^ ^ Mjup for the secondary, near the 
deuterium-burning mass limit. At the distance of 12.7 pc 
VHS 1256-1257 b is among the nearest currently known 
planetary mass companions detected by direct imaging. 
Moreover, it is one of the very few young, extremely red 
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L dwarfs with age constrained within a narrow range, 
given by the likely belonging to the Local Association 
and the absence of Li i in the primary. 

The Teff of ~ 900 K determined from evolutionary 
models based on the luminosity does not seem to be con¬ 
sistent with the expected Tes range of field dwarfs of 
similar spectral type and with the absence of methane, 
which is expected to appear in atmosphere cooler than 
1400 K. Following 
tion of clouds wit 


Barman et al. (2011a|b) the forma- 
ti substantial vertical thickness and 
non-equilibrium chemistry in a low-gravity object like 
VHS 1256-1257 b could provide an explanation of the 
apparent high atmospheric temperature (>1500 K) as 
compared to cooling track effective temperature predic¬ 
tions (900-1000 K). As atmospheric clouds are composed 
mostly of Fe and Mg-Si grains, we argue that the infall of 
planetesimals into this, and other young planetary mass 
objects, may enrich their atmospheres with key ingredi¬ 
ents for substantial cloud formation. As time progresses, 
the amount of infalling material will decrease, the sur¬ 
face gravity of the planet will increase, and therefore the 
formation of thick clouds will become less important. As 
a consequence, the majority of field mid L dwarfs will 
not display the extreme IR colors found in these young 
objects. An alternative hypothesis to aid the explanation 
of the very red colors involves presence of a warm debris 
disk or a dust-shell surrounding the object and causing 
the extinction. 

Since it is a relatively nearby and bright object near 
the deuterium-burning limit, it becomes one of the most 
promising targets to study the application of the deu¬ 
terium test. From the masses and separation of the 
components we estimate the orbital period to be about 
3900 yr. Assuming a circular orbit with a face-on ori¬ 
entation, the displacement caused by the orbital motion 
would be from 4 to 13 mas/yr, which will become feasible 
to measure in the next few years using precise astromet¬ 
ric observations. 
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